On the GNSS-R Interferometric Complex Field Coherence Time 
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O, 

In this paper we focus on the microwave bistatic scattering process, with the aim of deriving an 
expression for the interferometric complex field auto-correlation function from a static platform. 
We start from the Fresnel integral and derive the auto-correlation function in the Fraunhofer and 
Modified Fraunhofer regime. The autocorrelation function at short times can be expressed as a 
Gaussian with a direction dependent time scale. The directional modulation is a function of the 
angle between the scattering direction and the wave direction. The obtained relation can be used 
for directional sea state estimation using one or more GNSS-R coastal receivers. 
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1. Introduction 



The use of Global Navigation Satellite System (GNSS) signals reflected by the sea surface as 
a remote sensing tool has generated considerable attention for over a decade. Among several ap- 
plications, two classes have rapidly emerged in the community: sea-surface altimetry, which aims 
at retrieving the mean sea level like classical radar altimeters do, and sea-surface scatterometry 
or "speculometry" (a better term, as the GNSS-R return is dominated by quasi-specular scatter- 
ing) for the determination of sea roughness and near surface wind. This paper addresses a new 
application: sea state monitoring from low altitude, static platforms. 

The Oceanpal project at Starlab focuses on the development of technologies for operational 
in-situ or low-altitude water surface monitoring using GNSS-R. Oceanpal^ is a "down-to-Earth" 
offspring of technology developed within several ESA/ESTEC projects targeted on the exploita- 
tion of GNSS Reflections from space, following the proposal of Martin-Neira in 1993 |6|. This in- 
strument is to provide low cost/low maintenance sea state and altimetry measurements for coastal 
applications with the inherent precision of GNSS technology. Altimetric applications from the 
ground and from space have already been demonstrated in other campaigns . as well as 
airborne sea-state applications |5). 

The present paper was motivated by a recent Oceanpal experiment conducted by Starlab in 
collaboration with the Barcelona Port Authority for the demonstration of GNSS-R sea state moni- 
toring from low altitude, static platforms [9|. The HOPE 2 (Harbor Oceanpal Experiment) experi- 
mental campaign gathered GNSS-R data from a ground platform in a range of sea-state conditions. 
It was carried out during December 2003 at the Meteorological station of Porta Coeli belonging 
to the Barcelona Port Authority, located on the port breakers. Two antennas were deployed at ap- 
proximately 25 meters over the open sea surface to gather GPS-L1 signals. As usual in GNSS-R, 
one antenna was zenith-looking (Right Hand Circularly Polarized) to get the direct signal, and the 
other was oriented towards the open sea surface with the complementary polarization (Left Hand 
Circularly Polarization) to get the reflected signals. The antennas were connected to the instru- 
ment Front End, which then generated a digital bit-stream of unprocessed GPS data at IF. The IF 
data were recorded at a sample frequency of ~16 MHz, processed on board to produce complex 
waveforms, and transferred to the instrument mass storage unit. 

Our analysis for sea state begins with the interferometric complex field (ICF), defined at time 
t by Fj(t) = Ffi(t)/Frj(t), where Frj and Fr are the complex values at the amplitude peaks of 
the direct and reflected complex waveforms, respectively. The direct signal is thus used as a ref- 
erence to remove features unrelated to ocean motion, such as any residual Doppler, the navigation 
bit phase offset, or direct signal power variability. From a static platform, the ICF time history 
contains very valuable information on the sea state. 

The goal of the present analysis is to relate dynamics of the ICF to sea-surface geophysical 
parameters. Given the static character of the experiment and the very small scattering area — 
we recall the instrument was deployed at low altitude (~ 25 m) in this experiment — sea-surface 
roughness parameters such as the Directional Mean Square Slope (DMSS) are not foreseen to 
be quantifiable through the estimation of the width and orientation of the scattering area (Delay- 
Doppler mapping) as they were in — especially given the coarse nature of the C/A code which 
is available today. 

In this note we formulate an expression for the ICF autocorrelation function in the Fraunhofer 

' Patent pending. 
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and modified Fraunhofer approximations starting from the Fresnel integral for the scattered field 
(scalar physical optics). 

2. Fresnel integral 

We start our derivation from the Fresnel integral approximation for the scattered field and 
ignore for the time being curvature related polarization effects (see p. 380, Eq. 17 in |2|). The 
context is that of GPS reflections in a static receiver scenario. Although the GPS signal is not 
monochromatic, it can be thus approximated, since the code modulation (here represented by a[t\) 
adds relatively small bandwidth to the signal (co m << uj). This is an important consideration, 
since it means that there is no need to work with wave-packets. 

We start from a signal of the form a[t] exp[iuit] exp[m7r], as emitted by the satellite (the last 
term is the navigation bit, and we ignore Doppler and emission gain factors). The direct electric 
field at the surface is given by (refer to Figure 

—iujt+imr iks 

£?/(*) = a[t - s/c] , (2.1) 

4ns 

and the reflected electric field measured at the receiver is given by 

• _ '/ h' { 7' —J— S ) 

E(t) = ^ / VQTZa[t - (r + s) /c] e - iu > t+i ™- (q ■ h) dS, (2.2) 

An J rs 

where Q is the antenna gain, uj is the carrier frequency (including Doppler offsets), nn the (GPS) 
navigation bit, 1Z is the Fresnel coefficient, k = 2ir/\, with A «19 cm in GPS LI, r (s) is the 
distance between the receiver (transmitter) and each point of the sea-surface, h the normal to 
the surface, and q = (<fj_, q z ) is the scattering vector (the vector normal to the plane that would 
specularly reflect the wave in the receiver direction). This vector is a function of the incoming and 
outgoing unit vectors hi and h s , q = k(hi — h s ). 

Equation 12.21 iust states that each point in the surface acts like a local mirror source of the 
incoming field, incorporating the appropriate delay: the resulting field is the superposition of all 
those fields modulated by the antenna gain. Correlation of this signal with a replica of the "clean" 
direct signal will further modulate the support of this superposition, as we now discuss. 

After correlation of the field signal -f with an optimized code-carrier replica (such as a[t — 

t] exp(— i(u + 5uj)t]) we obtain [ 10 1, 



,„ -iAujt+inir p p ik(r+s) 



Fit) = / MTZ (q ■ h) dS, (2.3) 

4-7T J rs 



where Auj is the residual carrier frequency. The term M. = j\A(p,z) = \JQ{p, z) x(f>, z) rep- 
resents now here the antenna gain and Woodward Ambiguity Function (WAF) filtering on the 
surface. Here the WAF is given by 

1 pt+Ti 

X (p, t, 5t, 6u) = - a(t') a{t' + <5r) e -^ 4 ' (2.4) 

Jj Jt 



* Which at this stage has been transformed into a current by the antenna, then again to a voltage, down-converted, 
digitized, etc., facts that we ignore here as they are not crucial to the discussion. 



Starlab Technical Brief TB0005: ICF Coherence Time 3 

and 

St = t - (r + s)/c. (2.5) 

In the present case, \ can be approximated by [10] 

X (j?,t,8T,8f)=A(5T)-S(6u) (2.6) 

where 

l r T i 



l r 1 * 

A{5t) = — I a(t + t') a(t + t' + St)(U' 
m Jo 



1 -St/t c , \St\ < T c 

-r/Ti, \St\ > t c 



(2.7) 



and 

Sm . I T e^'cft' = ^y e-^^. (2.8) 

Note that this function has support near Slu rs 0, i.e., near \5f\ < l/2Tj. Similarly, A(<5r) provides 
most support to the integration area satisfying \t — (r + s)/c)\ < t c . Thus, F(t) effectively sums 
the field scattered from the surface at locations supported by the mentioned delay and Doppler 
conditions. The WAF's role is to simply modulate the support of the Fresnel integral. 

The GNSS-R receiver measures the (real part) of the electric field (including the carrier) 
through the coupling of the electric field to the antenna, which induces currents and voltage vari- 
ations. This voltage is a direct measure of the electric field (polarization issues aside), up to a 
multiplicative constant, which we ignore here — and noise. GNSS signal processing (the corre- 
lation process) allows for removal of carrier and code modulation and thus the recovery of the 
electric field spatial and slow temporal amplitude and phase variation as filtered by the WAF (as 
is nicely described in fTOll ). This is what Equation 12. 3l represents. 

We assume here that q ■ n ss k (small slope approximation and a small patch approxima- 
tion, with scattering and/or support only near the geometric flat surface specular point, or high- 
frequency limit). We first note that§ (Figure Q 

r' + s' = -q. (2.9) 

We now write s = s' + Sr and, for the emitter far field case which applies here (the emitting GPS 
satellite is very far in comparison with the scattering region size), 



s 




s' 2 + (Sr) 2 + 2s' • 5y w s'\ll + 2^-—^ w s' + s' ■ Sr. (2.10) 



Hence, we can rewrite Equation 12. 31 as 



■ u —iAuit+imr iks' r ik(r+s' -Sr) 

F(t)^—^— A e \ MR- dS. (2.11) 

47T s' J r 



*Asa convention, general vectors will be in bold face, while vectors laying in the plane will be displayed with an 
arrow on top. A hat will always indicate a unit vector. We will also write, for any vector, w = w± + w z z. 
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Transmitter 




Figure 1 : Definition of vectors used in analysis. The specular point is located at the origin. 



Now, the incoming field at the receiver is given by 

I(t) 



„— iAuit+iniv p ik\\s' — r'|| p iks' — ks'-r' 



Airs' 



Airs' 



(2.12) 



since | |s — r 1 1 « s — § • r . We can now write the emitter far field expression for the ICF (F/I) as 



Fi(t) « -ike lks r / MR 
and writing r' « r — Hz/r , we finally have 

F T (t) « -iifee 



ifc(r+s'-<5r) 



■d5, 



(2.13) 



p ik(r — Hz/r +s'-Sr) 
iks-r I Mn ^ 



(2.14) 



The exponent has now been expanded to first order in z, and the Gaussian approximation for the 
correlation function can now be used. We can further approximate this result by writing H/r m 
sin e and s' • 5r ~ s^_ • p — sin e 



F / (t)«-ifce lfcsr MR 

We now expand the r term in powers of Ap, 

r = \\r' + p\\ 2 



o ik(r — 2zsine+s' ± -p) 



dS. 



(2.15) 



(2.16) 



2r' 



r' + fW+T^'-L -p)' Z . 
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Keeping only first order terms leads to the so-called Fraunhofer approximation. A step beyond the 
Fraunhofer approximation is obtained by keeping the second order term. We call this the Modified 
Fraunhofer approximation. 

Keeping the second order term we can then write 

FAt) « -ik e% S r ^ />^(-2^sine-g x .p+^7(f' ± p) 2 ) d5 . (2 .1 7) 

r J 

and finally (§' • r' + r' = r' ■ (§' + f ) = 2iH sin e) 

p i2kH sine /■ to 

F/(t) « -ik- / Mne i(-^zsine-q ± -p+ — (r' ± -p) ) ^ (2 _ lg) 

3. Fraunhofer Approximation 

As an exercise, we carry out first the Fraunhofer approximation to the Fresnel integrand. We 
will see later that this approximation is not capable of representing directional dependence. 
Dropping the second order term in Equation l2.18l we get 



■i2kH sin e 



p iznn sm e /> 

Fi(t) = -ik- / MKe-^^dS. 



(3.1) 



Now we can compute the ICF autocorrelation function, 

T(At) = (dL f MM'KK' e~ i(l - Ai5r) dSdS'}, (3.2) 

and writing Sr = [p, z], and q = [q±, q z ] 

-k 2 r 

r(At) = -J- / MM'nW ( e -^-^p-i^ dS ds', (3.3) 



which (using q z = 2k sin e) we can rewrite — assuming Gaussian statistics for the surface to write 

ID 

/ -2iksine[z(p',t)-z(p',t+At)]\ _ -ik 2 sin 2 ea 2 [l-C(Ap,At)] to a\ 



where <j z is the standard deviation of the surface elevation and C(Ap, At) the spatio-temporal 
autocorrelation function of the surface — as 

7 9 

r(Ai) = —^- f MM'Wl' e ^x-Ap-4fc 2 sin 2 e^[ 1 _ C(A p-At)] dSd gl (35) 

We ignore, for the time being, WAF, antenna pattern and Fresnel coefficients space dependence to 
write 

r(At) « ^- f e ^ A ^ 4fc2 ^eal[i-c{ApAt)\ dSdS >_ (3 6) 

Ignoring these terms is licit at small At: then the integrand is strongly suppressed by the autocor- 
relation function itself. Nonetheless, in the next section we will keep them for completeness and 
clarity. 
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Now we write a simple Gaussian expression for the autocorrelation function and use the real 
version of the stationary phase approximation, the saddle-point approximation, to keep terms in 
A/9 only to second order, 



C(Ap, At) = e 2t ; 



--ApMAp 



e *% (1 - ApMAp) , 



(3.7) 



where 



~2 R ^' 


l// c 2 




R-ijj 






1 


cos ip sin ip 




\yii 







2 


— <i 


dn ip cos ip 







i/ii 





(3.8) 



cos ip — sin ip 
sin ^ cos ip 



and where /„ and l c are the up and cross "wave" correlation lengths, and ip the direction associated 
to /„. Another way to write M (using Dirac notation) is 



M = ^\v)(v\+^\v ± )(v ± \. 



(3.9) 



This expression for the autocorrelation function is directional, but does not include a surface 
current — a key aspect to see directional effects in the ICF, as we will see later. Now, define 



Then we use 



to calculate 



E[At] =4k z sm z eaie 2r i . 



-xAx+bxm _ -bA^ 1 b/i 



7f" 



detA' 



(3.10) 
(3.11) 



r(At) 



ATL 2 e 



„/2 



e iq ± -Ap+4k 2 sin 2 eai[C(Ap,At)] dSdS' 



AK 2 e- 



M 



,iq ± -Ap+S[At][l-ApMAp\ ^5^5" 



AH 2 e~ 4k sin e °z +H [ At l 



„/2 



iq ± .Ap-B.[At]ApMAp^g^gl 



An Z e -±k* S ir?ecrl+~{At](l-{q 1 _)M-lqj2) 



A ^ n 2 e -ik 2 S in 2 e<Tl+B.[At\(l-(q A _)M- l qj2) 

\E[At]\r' 2 



7T Z 



det (HfAijM- 1 ) 



det(M). 



(3.12) 
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Now, many approximations have been used so far. The important thing to retain is that the 
autocorrelation function is suppressed by the SWH dependent term and that the coherence function 
can be written in the form, 

r(At) « Ae -^sm^^ e -3[At](l-(g x )M-^ ± /2)_ (3J3) 

That is, suppressed by a SWH term, and with a temporal behavior modulated by the relationship 
of two directions: the surface autocorrelation direction, and the scattering axis. We can again look 
at short times and write and updated version of the autocorrelation function 

r(At)«ie ^F l Kq±) q±,) . (3.14) 

Since q± = in the near specular Fraunhofer scenario, we do not have a directional modification. 
This is the expression used in [9|. 

Directional sensitivity is present here only in the case the instrument is pointing away from the 
specular direction. In the specular scenario (e.g., the antenna is pointed to the geometric specular 
point, or the antenna gain is very low), directional sensitivity of the ICF coherence time will begin 
to appear only in the Modified Fraunhofer approximation and through the inclusion of a surface 
current term in the autocorrelation function, as we discuss next. 

4. Modified Fraunhofer 

A moment's thought will show that the only way that a modulating term associated to the 
scattering axis can appear is through the higher order terms in the Fresnel integral — the Fraunhofer 
limit is not sufficient. The modified Fraunhofer provides us with the reference direction r'± (the 
scattering axis). 

We begin the discussion starting with Equation 12. 181 

2ikH sine /■ 
Fl ( t ) a -ik- / Mn j(-2ksinez-q ± .p+^r(? x .p) ) ^ 

We will us a Gaussian term to simulate the impact of antenna gain and the WAF and write 

p 2ikH sine r , 

Fr(t) « -ikK- } / e i (- 2k ****-?±-ft-»p<.*±-i f ) +i P R P)dS, 



(4.1) 



where R is a regulating term to account for first chip zone, antenna gain, q ■ h, etc. We will write 

O^lrlXf'J, (4-2) 

so 

N = 7T1 I f ±)( f ±l + iR = + iR, (4.3) 

and 

2ikHs'me <• 

F^t) w -ikU- ; / e i( -~ 2k ^z-tx-P+fWp) dS. (4.4) 

We can now rewrite 

T(At) FS A ^ [ e^x-Ap-^fc 2 sin 2 ea'j[l-C(Ap ',At)]+ipNp-ip>N* p> ^SdS'. (4.5) 

r' 2 J 
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The next step is to define a wave direction. This is in fact defined already by the matrix M, 
but this directionality is not yet time-associated (which we need if we are to see the impact on the 
coherence time). We have to redefine the surface autocorrelation function to include the impact of 
something like a current. Now we write, 



C(Ap, At) = exp [- (Ap - vAt) M (Ap- vAt)] exp -A(Atj 

,2 



exp 



exp 



2/ 2 



-ApMAp- (At) 2 —, + Atv- Ap/li 



exp 



-A{At)' 



-At 2 ^ 
111 



1 - ApMAp + Atv- Ap/ll) exp -A(At) 



where we have used Mv = v/(2l 2 ). 

To be consistent with the earlier discussion, we rewrite this as 

-^£ f Atv ■ An 
C(Ap, At) = e ^F 1 - ApMAp + (3—^ 

V luT z 



(4.6) 



= e **?(!- ApMAp + p- Ap) 



(4.7) 



with p = [3 Atv / (I u t z ) . We expect << 1. Comparing this equation with Equation 13.71 we see 
we now have a new term, the spatio-temporal cross term. This expression for the autocorrelation 
function says that if an observer moves along in the "current" reference frame the surface will 
appear to decorrelate more slowly. 
We now have 



r(At) 



atz 2 

ATI 2 

r >2 ' 

An 2 



-ik 



.,. s'm 2 ea 2 / iq ± -Ap+4k 2 sin 2 e C (Ap,At)+ipN p-ip' N* p' jn jrf 



4k 2 sin 2 eu 2 / icf ± .Ap~+B[At](l-ApMAfa-p-Ap)+ip~Np*-ip'N*p' jhjc' 



4fc 2 sin 2 eo- 2 +S[At] 



Define now 



hence 



/' 



(f±-i£[At]p)-Ap-E[At](Ap~MAp)+ipNp*-if?N*p' ^SdS' 



Pffi = P + P, 
Pe =P~ P, 



Now, 






pNp- p'N*p' = - {p 9 (iR)p e + p e (iR)p e + P^Opq) 



(4.8) 
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A1Z _AU2 t! , ri 2 c „2 icta+I f lift, _in.rAtlflU.2f3 



r(At) ~ c -4fc 2 sin 2 £ g;+5[At1 /" g J(g ± -i_[At]p+-f^)p e -±p e (2E[At]M+R)p e -±PgjBp e ^ ^ 

We now use Equation l3.HI We compute first the "+" integral (with b = iOp e /2 and A = R/2), 
and then the "-" one, 



V(\f\ ~ ^_L-4fc 2 sin 2 e a 2 + ~[At] f i(q ± -i.E[At]p)P e -§P e (2S[At]A/+R+0i?- 1 0)p e ./____!__ JO 
1 j r' 2 7 " V det ( fi / 2 ) 

HAt) « 2?rA^ 2 p _ 4 fc 2 s in 2 e t 7 2 +5[Af] /" i(g ± -iH[At]p)p Q -ip e (2S[At]M+iJ+0iJ- 1 0)p e dfi , 

V ; ~ r'Vdet(i?) 7 e ' 

(4.9) 

and, finally 

rCAil ~ ^^^ p -4k 2 sm 2 ea 2 z +E[At]-^q ± -iS[At]p)-(2E[At]M+R+OR- 1 Oy 1 iq ± ^iE[At]p) 

1 j ~ r' 2 7ditCR) 



47T 2 



det (23[Ai]M + -R + Oi?- 1 ©) 



_ 2ttATZ 2 e - 4 fc 2 S in 2 e ( rg+H[At]-|(gl-i3[At]p).Q- 1 -(gl-i5[At]p) /_____ 
r /2 v / det( J R) V det< 3 

(4.10) 
with Q = 2E[Ai]M + -R + OR^ 1 0. Finally, using <fj_ = for the specular situation, we have 

FfAi) « 27T^^ 2 4fc 2 sin 2 eg 2 +E[At]+ 1 B[Af frEQ-l. g / JItT^ 



4.1. Isotropic antenna gain 

Recall that 



o = S-l«i><-l. 



and write 

R = I/d 2 , (4.12) 

where d is a length scale set by the antenna gain. Let us write 

|fj_) = cose|n<^). (4.13) 

We can now compute the determinant and inverse of the 2d matrix 

Q = 2Z{At]M + R + OR- 1 0. (4.14) 

Since (f j_|f j_) = cos 2 e, 



2r' 



2 =R^I \u<t,)(u<f,\ (4-15) 
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10 



and, also using the definition of M, 



[At] r . r . E[Ai] . ± ± I , dkcos 2 e \ ,..,., 



(4.16) 



S[At] IN,-/., /-[At] 1\ ,.!>,.,, (dkco$ 2 e 



+ Z?2 l«>H + 



+ ?l^l + 



^)(^| 



1.2 d 2 / ' /N V Z 2 ■ ,/2y ' ' - ! y 2 /-' 
Now we expand | u^) =\ v)(v \ u&)+ \ v )(v \ U&) = cos (p \ v) + simp \ v ) to write 

| iifj,) (u^ | = cos (p \ v)(v \ + sin <p | v ) (0 | + cos ip sin y> ( | v) (v | + | v }(v \) . (4. 17) 

We now define the scale 

2r' 
k = -n 5- ( 4 - 18 ) 



dk cos 2 e 

In the case of interest (a coastal application using the GPS C/A code in non-grazing angles, l v « 
1. Hence 

/-[At] 1 cos 2 cp\,*w^, /HfAt] 1 sin 2 w\ , ... , . 



cos 09 sm 09 / „ | ^ | \ 

^ — 72 — ^ (I v >< u \ + \ v ><*!)■ 



(4.19) 



Note that y> is the angle between the scattering and wave direction. It is straightforward to write 
the inverse of this matrix, using 



a b 
c d 



1 



ad — cb 



d -b 
—c a 



as 



det Q ■ Q 



-i (E[Ai\ 1 sin 2 09\ |AWA| (E\At\ 1 cos 2 ip 



cos ip sin p 
P 



(ivX^I + I^Xvl). 



For large d we ignore the middle term. 



detQ«^ + ^(^ + C ° sV 



Z 2 / 2 

S[At] 
Z 2 Z 2 ' 



/ 2 



/ 2 



n 



(4.20) 



u-'-Xu- 



(4.21) 



(4.22) 



In the last step we assume a small degree of spatial anisotropy in the ocean spectrum (/ = l u « l c ) 
and we use /^ << 1. 
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Finally 



Q 



-i 



2 I 2 sin^ ip \ . . / 2 , l 2 cos 2 ip' 



K + —< 



v){v\ + \li + — 



l2cOS ^ [nV (\v)(^\ + \^)(v\). 

f I 2 COS 2 </A aJ _ 



tr 1 }^ 



•[At] 



Z 2 sin 2 if 



v)(v 



•[At] 



v^iv- 



l 2 cos </? sin </? 



(|t5)(t)- L | + |fr L >(0|)- 



(4.23) 



Finally, returning to Equation l4.11l we can write 

2ttAK 2 



r(At) 



4fc 2 sin 2 e a2+H[At] + i^ = [At] 2 ^-5iS^-^ 



r' 2 v /det( J R) 



r z ' H|At] ' T Z 



/ 4vr 2 
detQ 



2TTATZ 2 d 2 lJ -4fc 2 sin 2 e CT 2 +E[At](l+ 



Z (Aty . 2 



„/2 



/ 4vr 2 
PI 



27T ATZ 2 d 2 l v l -4fc 2 sin 2 e ff 2 +4fc 2 sin 2 e CT 2 e ^Jfl+^i^ 

- e V T z 



j-sin <p 



„/2 



4vr 2 



\ 4fc 2 



sin 2 ecr 2 e 2t * 



4lT 2 AK 2 d 2 lJ -4fc 2 sin 2 e <7 2 +4A ; 2 sm 2 e <7 2 e " 2T Ffl+^^si 
- R V r z 



r' 2 2/csineo". 



Sln 2 ^ Ml 



4TT 2 AK 2 d 2 lJ -4fc 2 S in 2 eCT 2 [l-(l-|j)(l+f ^fsinV) 



r' 2 2/csme<T^ 

4yr 2 A^ 2 d 2 / ! ,/ -4fc 2 

c 

r /2 2fcsineov 



^(l-^sinV) 



(4.24) 



Now we can see that the combination of the Modified Fraunhofer expression for the field, and the 
new autocorrelation function bring two new directions to play: the scattering axis, and the wave 
direction. To first order, and allowing for empirical adjustment to our approximations, the new 
expected result will be 

^■sm 2 t(l-f3 2 sm 2 (p) 



T(At)^Ae 



(4.25) 



The coherence time of the ICF is now given by the width (second order moment) of this Gaussian 
function, 



A 



TF 



2ka z sin e J\ - (3 2 sin 2 tp vr sin e t/l - (3 2 sin 2 <p S WH 



(4.26) 



According to this model, tf depends on the electromagnetic wavelength and the ratio between 
the correlation time of the surface and the significant wave height (an inverse velocity). It is 
also apparent that the ICF is more coherent when the scattering direction and wave direction are 
perpendicular. 
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Equation 14.261 should provide the basis to determine direction in the waves (Oceanpal S WH 
Algorithm 2). Using data from visible satellites at any given time (with different elevations and 
azimuths) the optimal wave direction <p u and ocean "z-velocity" Z v = SWH/r^ parameters can be 
searched for (as described, e.g., in 0). We first write ip = <j) — 4> v , i.e., as the difference between 
the satellite azimuth and the wave direction. The we write, for each satellite in view, 



4 = = . ^r • ( 4 - 27 ) 



tt sin e k ^1- f3 2 sin 2 (<p k -<j) u ) Z " 

To make the measurement more robust, several Oceanpal instruments can be deployed in the area 
of interest. We then write and equation for each satellite-receiver link (Oceanpal SWH Algorithm 

2): 

rf = . A = i-. (4.28) 

vr sin e l ' k Jl - (3 2 sm 2 (4> % > k - (j) u ) ^ 

Finally, recall that in 0, based on the Elfouhaily et al. spectrum we derived a linear relation- 
ship between t z and the SWH: t z = a s + b s * SWH (with a s =0.167, 6 S =0.388, and an error of 
0.03 s). This relation turns out to be rather independent of wave age. Using it, we can now rewrite 
Equation l4.28l as 

j,k _ A a s + 6 s *SWH 



vr sin e*> fc Jl - (3 2 sin 2 (^. fc - <f> u ) SWH 



5. Conclusion 

In this note we have derived an expression which can be used for the definition of a semi- 
empirical algorithms for low altitude, static GNSS-R sea state and wave direction (as in [9|). 
The expression relates the interferometric field autocorrelation function coherence time with sea 
state parameters. Although directional sensitivity is already present in the Fraunhofer regime 
analysis with a moderate gain antenna factor pointing away from the specular, surface current 
sensitivity appears only in the Modified Fraunhofer regime. The derived expression is sensitive 
to the wave direction relative to the scattering direction and to the ocean "z-velocity", the ratio 
of SWH to surface coherence time, and can be used for coastal sea state monitoring using one or 
more receivers. 
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